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We derive a general expression for the cumulant generating function (CGF) of non-Markovian 
quantum stochastic transport processes. The long-time limit of the CGF is determined by a single 
dominating pole of the resolvent of the memory kernel from which we extract the zero- frequency 
cumulants of the current using a recursive scheme. The finite-frequency noise is expressed not only in 
terms of the resolvent, but also initial system-environment correlations. As an illustrative example 
we consider electron transport through a dissipative double quantum dot for which we study the 
effects of dissipation on the zero-frequency cumulants of high orders and the finite-frequency noise. 
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Full counting statistics (FCS) has recently attracted 
intensive theoretical [l[ and experimental Q attention. 
The interest stems from the usefulness of FCS as a 
sensitive diagnostic tool of stochastic electron transport 
through mesoscopic systems. Detectable mechanisms in- 
clude quantum- mechanical coherence, entanglement, dis- 
order, and dissipation Mathematically, FCS en- 
codes the complete knowledge of the probability distri- 
bution Pin, t) of the transmitted number n of electrons 
or, equivalently, of all corresponding cumulants. Non- 
zero higher order cumulants describe non-Gaussian be- 
havior. The study of counting statistics for stochastic 
processes in general is of broad relevance for a wide class 
of problems, also outside mesoscopic physics. For exam- 
ple, rare events, whose study has become an important 
topic within non-equilibrium statistics of stochastic sys- 
tems in physics, chemistry, and biology are reflected 
in higher order cumulants. Efficient methods for evalu- 
ating the counting statistics of stochastic processes are 
therefore of urgent need. 

In this Letter, we present a method which unifies and 
extends a number of earlier approaches to FCS within a 
generalized master equation (GME) formulation 
The earlier approaches have in practice been limited to 
systems with only a few states 0, @, or only the first 
few current cumulants [5]. In contrast, our theory en- 
ables studies of a much larger class of problems: Eval- 
uation of zero-frequency current cumulants of very high 
orders for non-Markovian systems with many states is 
now possible. Furthermore, the method allows us to 
develop a general approach to the finite-frequency cur- 
rent noise of non-Markovian transport processes. In the 
case of finite-frequency noise, we show that not only the 
memory kernel but also initial system-environment cor- 
relations are crucial. Such correlations can be, and have 
been 0] , neglected for non-Markovian processes at low 



frequencies, but must be included at frequencies compa- 
rable with the internal frequencies of the system. We 
demonstrate our methods on a system of recent experi- 
mental relevance, namely transport through a dissipative 
double quantum dot @, H| , but they may easily be ap- 
plied to other electronic (or photonic) counting systems, 
as well as other counted quantities, such as heat or work, 
in non-equilibrium systems Q. 

Non-Markovian GME. Consider a nanoscale transport 
system governed by a generic non-Markovian GME of the 

form da, El 



(1) 

Here, the reduced density matrix of the system p{t) 
has been resolved into components p(n, t) correspond- 
ing to the number of electrons n passing through the 
nanosystem within time-span [0, i]. The memory ker- 
nel W describes the influence of the environment on the 
dynamics of the system, while the inhomogeneity 7 ac- 
counts for initial correlations between system and en- 
vironment. Both VV and 7 decay with time, usually 
on a comparable timescale, so that 7 is irrelevant for 
the long-time limit. The inhomogeneity 7 does, how- 
ever, play a crucial role at finite times. The probabil- 
ity distribution for the number of transferred charges is 
P(n,t) = Tr{p(n, t)\. The corresponding cumulant gen- 
erating function (CGF) S(\,t) is defined as e 5 ^'*-* = 
J2 n P( n : t)e mx . In Laplace space Eq. |T]) leads to the al- 
gebraic expression p(x, z) = G(x,z)[p(x, t = O)+7(x,*0], 
where Q(x> z ) = [z— W(x, z)] -1 is the resolvent of the ker- 
nel, and p(x, z) = Yin Jo°° dtp(n,t)e mx ~ zt and similarly 
for 7(x7 z ) and W(x,z). Inverting the Laplace transfer- 
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mation, the CGF then becomes 



2ni 



] = TT- / dz{Q{ X ,z))e z \ 



(2) 



where a is a real number, chosen such that all singulari- 
ties of the integrand are situated to the left of the verti- 
cal line of integration. We have moreover introduced the 
notation (G( X ,z)) = Tr{£?(x, z)[p( X ,t = 0) + j(x,z)}}. 
Equation |(5J) contains the full statistical information 
about the charge transfer process. It is a powerful formal 
result, but it also leads to useful practical schemes, as we 
shall now demonstrate. 

Zero-frequency FCS. Consider first the zero-frequency 
cumulants of the current, defined as ((I m )} = 
Tt °G>(fx) X " ) Ix^o.t^oo, TO = 1,2, .... We assume that the 
system with the counting field X set to zero tends ex- 
ponentially to a unique stationary state determined by 
the 1/z pole of the resolvent Q(x = 0, z). The sta- 
tionary state is given by the eigenvector corresponding 
to the zero-eigenvalue of W = W( X = 0, z = 0), i.e., 
limt_>oo = p stat , where p stat is the normalized solu- 
tion to W/5 stat = 0. With finite values of x, an eigenvalue 
Ao(x, z) develops adiabatically from the zero-eigenvalue, 
such that Ao(0, z) = 0, and the long-time behavior is de- 
termined by the isolated pole structure l/[z — Xq(x, z)] 
of Q(x, z) close to zero. This pole zq( X ) solves 



zo ~ A (x, z ) = 0, 



(3) 



and goes to zero with x going to zero, i.e., zq(0) = 0. We 
thus find e si - x ^ ->■ V{ X )e Zo( - x '> t for large t, where V( X ) 
is a time-independent function depending on the initial 
conditions and correlations. The current cumulants then 
read ((i™)} = ix^o- In the Markovian limit for 

the kernel W( X , z — > 0) we get z {x) = Ao(x> 0) as found 
in Refs. [IS]. 

Recursive scheme. When the involved matrices are 
large, it may be non-trivial to determine the full X and 
z dependence of the eigenvalue X ( X ,z), and thereafter 
solve Eq. ([3|). Instead, we expand the eigenvalue as 

A ( x , z) = J2Z=o {j TTir c{k ' l) with c(0 ' L) = °. and calcu " 
late the expansion coefficients recursively using Rayleigh- 
Schrodinger perturbation theory 



l — n V / n V / 
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with K,L = 0,1,2,.... Here, ((0| solves ((0|W = 
0, while |0( 0,0 '}) is the stationary state p stat . More- 
over, W(x, z) = W(x, z) — W has been expanded as 
Wfc*) = E^o^ffW^ with W im = 0. Fi- 
nally, the pseudoinverse of the kernel is 1Z = QW 1 Q 



with Q = l- |0( ' '})((0| [21]. With the c^-^'s at hand 
we can solve Eq. ([3]) for z { X ) to a given order in x, an d 
from the expansion z (x) = SnLi {(I n )) extract the 
zero-frequency cumulants of the current: 

1 1 



((I N )) 

p{K,L) 



N 

m V ±Lp(N-k,i) jk,i) 

2^ kin 

k,l=0 



E i(I n )) p(K-n,L-l) 



(5) 



n=l 



with L = 0, 1, 2, . . ., and K,N =1,2,3,.. .. For the aux- 
iliary quantity p( K > L ), we have P^ K ^ = 5 K ,o, P (0 ' L) = 
Sox, and P( K ~V = 0. 

We illuminate the recursive scheme by evaluating the 
first three cumulants of the current using Eq. ([5]); the 
mean current, the variance (the noise), and the skewness: 

((I 1 )) =c (1 '°>, 

((I 2 )) = C ( 2 »°) + 2c< 1 '°>c< 1 ' 1 > ) 

«/ 3 )) - c ( 3 '°) + 3c( 2 < ( V 1 < 1 > (6) 
+ 3c(^°) UWcQA + 2( C ( 1 ' 1 )) 2 + c^' . 

Higher order cumulants can be obtained in a simi- 
lar manner, analytically or numerically. Coefficients 
of the form c^' ) are purely Markovian quantities, 
and the mean current is thus not sensitive to non- 
Markovian effects, whereas higher order cumulants are 
6, 14]. From Eq. (gj) we find for the coefficients c (K - L) , 
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= <<o|W (1 ' n) |o (0 ' 0) )>, = ((o|(W 



e.g.. 

W (1,0) ftW (0,1) )|o(°'°))), 
2W il ' 0) nw im 



(i,i) 



and = ((0|(W (2 ' 0) 



)|0(°'°))). Evaluation of the pseudoin- 
verse 1Z amounts to solving matrix equations which is 
feasible even with very large matrices [13j. Numerically, 
the recursive scheme is stable for very high orders of cu- 
mulants (> 20) as we have tested on simple models. 

Double quantum dot. We illustrate our method by con- 
sidering a model of charge transport through a Coulomb 
blockaded double quantum dot in a dissipative environ- 
ment Q. Maximally one additional electron is allowed 
on the double quantum dot. The Hamiltonian of the 
double quantum dot is Hs = §Sz + T c s Xl where the 
pseudo-spin operators are s z = \L)(L\ — \R)(R\ and s x = 
\L)(R\ + \R)(L\, respectively. The tunnel coupling be- 
tween the two quantums dot levels \L) and \ R) is T Cl while 
e is their detuning. The pseudo-spin system is tunnel- 
coupled to left (L) and right (i?) leads via the tunnel- 
Hamiltonian H T = Yjk a ,a=L,R^Yk a c\ a |0)(a| +h.c), with 
both leads described as non- interacting fermions, i.e., 
H a = J2k a £ k a cl a Ck a ,ct = L,R. Dissipation is provided 
by a reservoir of non-interacting bosons that couple to 
the s z component of the pseudo-spin. The Hamiltonian 
is then H = Hs + H T + H L + H R + H B + V B s z , where 
Hb = J2j tv^jd]d and V B = J2j + a j)- 
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To describe charge transport through the system we 
trace out the leads following Gurvitz and Prager 
leading to an equation of motion for the reduced den- 
sity matrix a = (a 00 ,a LL ,a RR ,a LR ,a RL ) T of the dou- 
ble dot and the bath of bosons. The elements &ij are 
still operators in the Hilbert space of the boson bath. 
Charges are assumed to enter the left quantum dot from 
the left lead and leave from the right quantum dot via 
the right lead with energy-independent rates r Q (e) = 
^Ylk l^fcJ 2 ^( e — e O — T Q , a = L,R. This approach 
is valid to all orders in the tunnel coupling T c under the 
assumption of a large bias across the system [15 1. 

Next, we consider the electronic occupation probabil- 
ities pi = Trs{i7ii}, i = 0,L,R, where Tt r is a trace 
over the bosonic degrees of freedom. A closed system of 
equations is obtained by assuming that the boson bath 
at any time is in local equilibrium corresponding to the 
given charge state: an ~ pi <g> &i(f3), i — L,R, where 
a L/R ([3) = e-^/Tr B {e-^}, ff<±> = H B ± V B , 
and (3 = 1/fcnT is the inverse temperature (see e.g. Sec. 
IV C in Ref. [HI)- This approximation is valid when the 
bath-assisted hopping rates Tg(z) (proportional to T c 2 , 
see below) are much smaller than Tl/ R . The memory 
kernel for this model, with p = (po, p^, p R ) T , then reads 



W( X ,z) 




r 



R 

(-)/ 



B (*) 

\z)- 



(J) 



Here, the counting held x nas been introduced in 
the off-diagonal element containing the rate T R — > 
T R e lx , corresponding to counting of the number 
of electrons that have been collected in the right 
lead. The bath-assisted hopping rates entering the 
kernel are r^ ±} (z) = T 2 [g^{z±) + g { -\z T )] with 
gM{z) = f °°dte- w ^-*\ W(t) = I~dwJ(u>){[\ - 
cos(arf)] coth (/3cj/2) -I- i sm(cut)} / oj 2 , and z± = z ± is + 
T R /2. The spectral function of the heat bath is J(u>) = 
. \gj\ 2 5(uj — u>j), and below we show results for Ohmic 
dissipation, Jci(uj) — 2awe~"^ c , when the rates can be 
evaluated either analytically (for (3uj c 3> 1) or numeri- 
cally. 

In Fig. Q] we show the first three cumulants of the cur- 
rent as functions of the level detuning e with different 
dissipations strengths a. As an illustrative example the 
15'th cumulant of the current, ((I 15 )), is also shown. As 
the dissipation strength is increased, a clear suppression 
of the coherent features (with a = 0) is seen. The in- 
creased coupling to the heat bath tends to localize the 
electron to one of the two quantum dots, thereby sup- 
pressing the effects of the coherent coupling between 
them. As a result a cross-over from coherent to sequential 
tunneling is observed with increasing a. For large a's, the 
sequential tunneling process between the two quantum 
dots constitutes a 'bottle-neck' and the cumulants con- 
sequently approach the Poisson limit ((/ m )} / ((I 1 )) = 1, 
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FIG. 1: (color online). Zero-frequency cumulants of the cur- 
rent as functions of the level detuning e with different dis- 
sipation strengths a. Parameters are V = Tl — T_r = 0.5, 
T a = 0.1, k B T = 0, loc = 500, and 2ixa = (full line), 0.2 
(full line with dots), 0.5 (dashed line), 1 (dashed line with 
dots). A large bias is applied across the system. 



m = 1,2,3,.... The typical behavior of cumulants is, 
however, factorial growth, i.e., | ((/'"}} | ~ c q m ml for some 
constants c, q > 0. In the coherent case (a — 0), this 
behavior is clearly seen for the 15'th cumulant, demon- 
strating its high sensitivity to dephasing and decoherence 
mechanisms. 

Finite-frequency noise. The expression for the CGF, 
Eq. © , allows us also to study the finite- frequency spec- 
trum of the second cumulant of the current, the (sym- 
metrized) current noise [22] ]. expressed by MacDonald's 
formula as S n (u) = to J °° (it sin (ujt)((I 2 ))(t) [II, [13, El , 
where ({I 2 ))(t) = |^g*#|^ . We then find 



Sii{uj) = 



i) 2 



2 d(i X y 



[(Q(X> z = + (u! -> -cj)] | x _> < 

(8) 

In order to evaluate this expression, we need to choose 
p(x, t = 0) appropriately and find the inhomogeneity 
7(x, z) as they enter the definition of {Q{\,z)). Fol- 
lowing Ruskov and Korotkov [lj| we assume that the 
system evolves from to = —00, such that the elec- 
tronic occupation probabilities at t = 0, where electron 
counting begins, have reached the stationary state, i.e., 
p(n, t = 0) = (5 r i,o/S stat - For this model the inhomogeneity 
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FIG. 2: (color online). Finite-frequency current noise spec- 
trum for different temperatures (left figure) and dissipation 
strengths (right figure). Left figure: a — 0.005, fcgT = (full 
line), 1 (full line with dots), 2 (dashed line), 5 (dashed line 
with dots). Right figure: k B T = 0, and a = 0.01 (full line), 
0.02 (full line with dots), 0.03 (dashed line), 0.05 (dashed line 
with dots). Other parameters are Tl = T_r = 0.01, T c = 3, 
e = 10, A = + (2T C ) 2 , and lu c 
applied across the system. 



500. A large bias is 



is independent of the counting field [12j |: 
m = W-W(x = 0,*)^ at 



dissipation strength, which simultaneously reduces the 
effect of the coherent coupling. For larger a's (~ 0.1) a 
change of line shape is observed (not shown) . 

In conclusion, we have presented a general theory for 
current fluctuations in non-Markovian quantum trans- 
port systems. Our methods allow us to calculate recur- 
sively zero-frequency current cumulants of very high or- 
ders, governed by a single dominating pole of the resol- 
vent of the memory kernel, as well as the finite-frequency 
current noise, which is given not only by the resolvent, 
but also initial correlations. As an illustrative example 
of our approach, we have considered transport through a 
dissipative double quantum dot for which we have studied 
the effects of dissipation and temperature on the current 
cumulants of very high orders and the finite-frequency 
current noise. 
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of the Czech Science Foundation. 
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We see that the effects of the initial correlations ac- 
counted for by j(z) vanish in the long-time limit. 
Moreover, since Wp stat = 0, we find (Q(x,z)) = 
Tr{C*(x, z)Q~ 1 (x = 0, z)p stlit }/z from which we can cal- 
culate the finite-frequency current noise. We note that 
only the proper inclusion of the inhomogeneity ensures 
a correct finite-time behavior, such as proper normaliza- 
tion of Tr{p(t)} = (Q(x = 0,i)) = Tr{p stat } = 1 at all 
times. 

Displacement currents (due to finite capacitances be- 
tween quantum dots and leads) can be included in the 
finite- freq uency current noise via the Ramo-Shockley the- 
orem 20]. Evaluating the full current noise spectrum 
then requires an additional counting field accounting for 
tunneling from the left lead [18| ■ For the shown results we 
have included displacement currents and assumed identi- 
cal capacitances between the left /right quantum dot and 
the left/right lead. In Fig. [5] we show the finite-frequency 
current noise Sh(uj) at frequencies around the hybridiza- 
tion energy A = e 2 + (2T C ) 2 , where signatures are ex- 
pected in the current noise spectrum 0, [13, Res- 
onances are not observed exactly at the 'bare' value 
uj = A, but are shifted towards lower frequencies. This 
renormalization occurs due to coupling to the heat bath 
which dresses the eigenspectrum of the electronic degrees 
of freedom, similar to the Lamb shift in atomic physics. 
The left figure shows how the signatures due to the coher- 
ent coupling between the two quantum dots are washed 
out with increasing temperature. In the right figure, we 
show how the frequency shift increases with increasing 
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